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ABSTRACT 

Context. For the spectral analysis of high-resolution and high-signal-to-noise (S/N) spectra of hot stars, state-of-the-art non-local 
thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic 
data that is used for their calculation. 

Aims. Reliable BaV- VII oscillator strengths are used to identify Ba lines in the spectra of the DA-type white dwarf G191-B2B and 
the DO-type white dwarf RE 0503-289 and to determine their photospheric Ba abundances. 

Methods. We newly calculated BaV- VII oscillator strengths to consider their radiative and collisional bound-bound transitions in 
detail in our NLTE stellar-atmosphere models for the analysis of Ba lines exhibited in high-resolution and high-S/N UV observations 
of G191-B2B and RE 0503-289. 

Results. For the first time, we identified highly ionized Ba in the spectra of hot white dwarfs. We detected Ba VI and Ba VII lines 
in the Far Ultraviolet Spectroscopic Explorer (FUSE) spectrum of RE 0503-289. The BaVl/BaVII ionization equilibrium is well 
reproduced with the previously determined effective temperature of 70 000 K and surface gravity of log g = 7.5. The Ba abundance is 
3.5 ± 0.5 x 10 -4 (mass fraction, about 23 000 times the solar value). In the FUSE spectrum of G191-B2B, we identified the strongest 
Ba VII line (at 993.41 A) only and determined a Ba abundance of 4.0 ± 0.5 x 10 -6 (about 265 times solar). 

Conclusions. Reliable measurements and calculations of atomic data are a pre-requisite for stellar-atmosphere modeling. Observed 
BaVI-VII line profiles in two white dwarfs’ (G191-B2B and RE 0503-289) far-ultraviolet spectra were well reproduced with our 
newly calculated oscillator strengths. This allowed to determine the photospheric Ba abundance of these two stars precisely. 

Key words, atomic data - line: identification - stars: abundances - stars: individual: G191-B2B - stars: individual: RE 0503-289 - 
virtual observatory tools 


1. Introduction 

In recent analyses of the hydrogen-rich DA-type white 
dwarf (WD) G191-B2B (effective temperature T e g = 60000, 
surface gravity log (g/cm/s 2 ) = 7.6, Rauch et al. 2012, 2013, 
2014) and the hydrogen-deficient DO-type WD RE 0503-289 
(7’gff = 70000, log g = 7.5, Werner et al. 2012; Rauch et al. 
2014), numerous lines of the trans-iron elements Zn, Ga, Ge, As, 
Se, Kr, Mo, Sn, Te, I, and Xe were identified. This substantially 
reduced the number of unidentified lines in the spectra of these 
two WDs. For precise abundances determinations, reliable tran- 
sition probabilities are mandatory - these are necessary not only 


* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy 
Inc., under NASA contract NAS5-26666. 

** Based on observations made with the NASA-CNES-CSA Far 
Ultraviolet Spectroscopic Explorer. 

*** Tables 1-3 are only available at the CDS via anonymous ftp to 
cdsarc .u-strasbg . fr (13® . 79 . 128 . 5) or via 
http : //cdsarc .u-strasbg . fr/viz-bin/qcat?J/A+A/566/Al® 


for the identified lines themselves but for the complete model 
atom that is considered in the model-atmosphere and spectral- 
energy-distribution (SED) calculations. Thus, abundance deter- 
minations were so far restricted to Zn (Rauch et al. 2014), Ge 
(Rauch et al. 2012), Kr, Xe (Werner et al. 2012), and Sn (Rauch 
et al. 2013). 

A close inspection of the still unidentified lines in the 
far-ultraviolet (FUV) spectrum of RE 0503-289 showed ab- 
sorption features at the locations of the strongest Bavi and 
BaVll lines as given by NIST 1 . Therefore, we calculated Bav- 
VII transition probabilities (Sect. 2) and employed our NLTE 2 
model-atmosphere code (Sect. 4) to perform test calculations 
(Sect. 5) to find the strongest Ba lines in the model. Then, 
we used these strategic lines to determine the Ba abundances 
of G191-B2B and RE 0503-289 and searched for other, weak 
Ba lines (Sect. 6). We summarize our results and conclude in 
Sect. 7. 


1 National Institute of Standards and Technology, http://www. 
nist . gov/pml/data/asd. cfm 

2 Non-local thermodynamic equilibrium. 
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2. Transition probabilities in BaV, BaVl, and BaVII 

Very few studies have been focused on the determination of 
electric dipole transition rates in BaV, BaVl, and BaVII so far. 
To our knowledge, the only available data were recently pub- 
lished by Sharma et al. (2014) and Sharma et al. (2013) for 
BaV and BaVII, respectively. More precisely, these authors re- 
ported oscillator strengths and transition probabilities computed 
using rather limited configuration interaction models based on 
the Hartree-Fock approach due to Cowan (1981) combined with 
a semi-empirical least-squares fit of radial energy parameters. In 
order to get a uniform set of oscillator strengths for all the tran- 
sitions of Ba ions observed in the present work, we decided to 
perform new calculations including a larger amount of electron 
correlation and hence improving the previous investigations of 
Sharma et al. (2013, 2014). The method adopted here was the 
relativistic Hartree-Fock (HFR) approach with core-polarization 
corrections (see e.g., Quinet et al. 1999, 2002). 

For BaV, configuration interaction was considered among 
the configurations 5s 2 5p 4 , 5s 2 5p 3 4f, 5s 2 5p 3 5f, 5s 2 5p 3 6f, 
5s 2 5p 3 7f, 5s 2 5p 3 6p, 5s 2 5p 3 7p, 5s 2 5p 2 4f 2 , 5s 2 5p 2 5d 2 , 5s 2 5p 2 6s 2 , 
5s 2 5p 2 6p 2 , 5s 2 5p 2 5d6s, 5s 2 5p 2 4f6p, 5s5p 4 5d, 5s5p 4 6d, 5s5p 4 6s, 
and 5p 6 for the even parity, and 5s5p 5 , 5s 2 5p 3 5d, 5s 2 5p 3 6d, 
5s 2 5p 3 7d, 5s 2 5p 3 6s, 5s 2 5p 3 7s, 5s 2 5p 2 4f5d, 5s 2 5p 2 4f6d, 

5s 2 5p 2 4f6s, 5s 2 5p 2 5d6p, 5s 2 5p 2 6s6p, 5s5p 4 4f, 5s5p 4 5f, 

5s5p 4 6f, and 5s5p 4 6p for the odd parity. The core-polarization 
parameters were the dipole polarizability of a BalX ionic core 
as reported by Fraga et al. (1976), i.e., a& = 0.54 a.u., and 
the cut-off radius corresponding to the HFR mean value ( r ) 
of the outermost core orbital (4d), i.e., r c = 0.79 a.u. Using 
experimental energy levels reported by Sharma et al. (2014), 
the radial integrals (average energy. Slater, spin-orbit and 
effective interaction parameters) of 5s 2 5p 4 , 5s 2 5p 3 6p, 5s5p 5 , 
5s 2 5p 3 5d, and 5s 2 5p 3 6s configurations were optimized by a 
well-established least-squares fitting procedure in which the 
mean deviations with experimental data were found to be equal 
to 185 cm -1 for the even parity and 217 cm -1 for the odd parity. 

For BaVl, the configurations included in the HFR model 
were 5s 2 5p 3 , 5s 2 5p 2 4f, 5s 2 5p 2 5f, 5s 2 5p 2 6f, 5s 2 5p 2 7f, 5s 2 5p 2 6p, 
5s 2 5p 2 7p, 5s 2 5p4f 2 , 5s 2 5p5d 2 , 5s 2 5p6s 2 , 5s 2 5p6p 2 , 5s 2 5p5d6s, 
5s 2 5p4f6p, 5s5p 3 5d, 5s5p 3 6d, 5s5p 3 6s, and 5p 5 for the odd 
parity, and 5s5p 4 , 5s 2 5p 2 5d, 5s 2 5p 2 6d, 5s 2 5p 2 7d, 5s 2 5p 2 6s, 
5s 2 5p 2 7s, 5s 2 5p4f5d, 5s 2 5p4f6d, 5s 2 5p4f6s, 5s 2 5p5d6p, 
5s 2 5p6s6p, 5s5p 3 4f, 5s5p 3 5f, 5s5p 3 6f, and 5s5p 3 6p for the even 
parity. In this ion, the semi-empirical process was performed 
to optimize the radial integrals corresponding to 5s 2 5p 3 , 5s5p 4 , 
5s 2 5p 2 5d, and 5s 2 5p 2 6s configurations using the experimental 
levels reported by Tauheed & Joshi (1994). The mean deviations 
between calculated and experimental energy levels were 5 cm -1 
and 128 cm -1 for odd and even parities, respectively. Core- 
polarization effects were estimated using the same ay and r c 
values as those considered in Ba V. 

Finally, a similar model was used in the case of BaVII for 
which the 5s 2 5p 2 , 5s 2 5p4f, 5s 2 5p5f, 5s 2 5p6f, 5s 2 5p7f, 5s 2 5p6p, 
5s 2 5p7p, 5s 2 4f 2 , 5s 2 5d 2 , 5s 2 6s 2 , 5s 2 6p 2 , 5s 2 5d6s, 5s 2 4f6p, 
5s5p 2 5d, 5s5p 2 6d, 5s5p 2 6s, and 5p 4 even-parity configurations 
and the 5s5p 3 , 5s 2 5p5d, 5s 2 5p6d, 5s 2 5p7d, 5s 2 5p6s, 5s 2 5p7s, 
5s 2 4f5d, 5s 2 4f6d, 5s 2 4f6s, 5s 2 5d6p, 5s 2 6s6p, 5s5p 2 4f, 5s5p 2 5f, 
5s5p 2 6f, and 5s5p 2 6p odd-parity configurations were explic- 
itly included in the HFR model. Here also, we used the same 
core-polarization parameters as those considered for BaV. The 
semi-empirical optimization process was carried out to adjust 
the radial parameters in 5s 2 5p 2 , 5p 4 , 5s 2 5p4f, 5s 2 5p6p, 5s5p 2 5d, 
5s5p 3 , 5s 2 5p 2 5d, and 5s 2 5p 2 6s with the experimental energy 


Table 4. Statistics of our Ba model atom. 


Ion 

NLTE levels 

LTE levels 

Lines 

Ba IV 

1 

6 

0 

V 

123 

0 

981 

VI 

47 

0 

162 

VII 

121 

0 

452 

VIII 

1 

0 

0 


levels classified by Sharma et al. (2013) below 350000 cm -1 and 
320000 cm -1 for even and odd parities, respectively. In fact, it 
was found that many levels above those limits overlap unknown 
levels and are strongly mixed with states belonging to higher 
configurations such as 5s 2 4f 2 (even parity), 5s5p 2 4f, 5s 2 4f5d, 
and 5s 2 5p6d (odd parity). It was then extremely difficult to es- 
tablish an unambiguous correspondence between the calculated 
values and the experimentally determined level energies. The use 
of radial parameters published by Sharma et al. (2013) was un- 
fortunately not of great help for making the identifications more 
reliable, the set of interacting configurations being not the same 
as the one considered in our work. For the levels considered in 
our fitting process, the mean deviations between calculated and 
experimental values were found to be equal to 140 cm -1 (even 
parity) and 196 cm -1 (odd parity). 

Figure 1 shows Grotrian diagrams of BaV- VII including all 
levels and transitions from Tables 1-3. 


3. Observations 

In this analysis, we use the FUSE 3 spectrum of RE 0503-289 
and the FUSE and HST/STIS 4 spectra of G 1 9 1 -B2B that are de- 
scribed in detail by Werner et al. (2012) and Rauch et al. (2013, 
2014), respectively. 

The FUSE spectrum covers the wavelength range 910 A < 
A < 1188 A. Its resolving power is R = A/AA ~ 20000. The 
HST/STIS spectrum of G191-B2B is co-added from 105 obser- 
vations with the highest resolution (grating E140H, R ~ 1 18 000, 
1 145 A < A < 1750 A) available via MAST. 

4. Model atmospheres and atomic data 

Our model atmospheres (plane-parallel, chemically homoge- 
neous, in radiative and hydrostatic equilibrium) were calculated 
with the Tubingen NLTE model-atmosphere package (TMAP 3 , 
Werner et al. 2003; Rauch & Deetjen 2003). 

The model atoms (including Ba, Table 4 displays its statis- 
tics) that are considered in our calculations are provided by the 
Tubingen Model-Atom Database (TMAD 6 , Rauch & Deetjen 
2003), that was constructed in the framework of the German 
Astrophysical Virtual Observatory (GAVO 7 ). The SEDs that 
were calculated for this analysis are available via the registered 
Theoretical Stellar Spectra Access (TheoSSA 8 ) VO service. 

3 Far Ultraviolet Spectroscopic Explorer. 

4 Hubble Space Telescope/Space Telescope Imaging Spectrograph, for 
our high-resolution spectrum of G191-B2B, see http://www. stsci . 
edu/hst/observatory/crds/calspec.html 

5 http : //astro . uni - tuebingen . de/~TMAP 

6 http : //astro . uni - tuebingen . de/~TMAD 

7 http://www.g-vo.org 

8 http://dc.g-vo.org/theossa 
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Fig. 1. Grotrian diagrams of our BaV (top), 
BaVI (middle), and BaVII (bottom) model 
ions. Horizontal bars indicate levels, gray 
lines represent radiative transitions with known 
/ values (Tables 1-3). Our strategic BaVI 
and BaVII lines in the FUSE wavelength 
range (Fig. 6) are labeled with their wave- 
lengths. The long-dashed, horizontal lines in 
each panel show the ionization energies of 
468 000 cm -1 , 573 000 cm -1 , and 694 000 cm -1 
of BaV, BaVI, and BaVII, respectively. 


The start models for our calculations are the most elaborated 
model atmospheres for both stars (Rauch et al. 2014). As our 
main model-atmosphere program would not compile if the ar- 
ray sizes were increased to accommodate the high number of 
atomic levels treated in NLTE and the high number of radiative 
and collisional transitions, (cf. Rauch et al. 2014), we simply 


reduced the number of the Zn levels treated in NLTE to one per 
ion to create a TMAP executable. Because the Zn opacities were 
already considered in detail in our start models, the atmospheric 
structure and the background opacities are well modeled. To cal- 
culate the NLTE occupation numbers of Ba, we performed line- 
formation calculations, i.e., we kept temperature and density 
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Fig. 2. Ba ionization fractions in our models for G191-B2B ( top panel, 
H+Ba) and RE 0503-289 ( bottom panel, He+Ba). m is the column 
mass, measured from the outer boundary of our model atmospheres. 


structure of our model-atmospheres fixed. Since the impact of 
Ba on the atmospheric structure of our H+Ba and He+Ba test 
models (Sect. 5) was found to be marginal, this is the best prac- 
tice. For the subsequent SED calculations, we considered the 
complete model ions and resumed the missing Zn occupation 
numbers from the start models. 


5. Preliminary analysis 



1000 1100 
wavelength / X 


Fig. 3. Synthetic spectra (normalized to the models’s continuum fluxes) 
of our H+Ba model for G191-B2B {top panel) and our He+Ba model 
for RE 0503-289 {bottom panel). All synthetic spectra displayed in this 
paper are convolved with Gaussians to simulate the instruments’ resolu- 
tion. Here, for FUSE, we used an FWHM (full width at half maximum) 

of 0.06 A. 


For a preliminary analysis, we calculated models that con- 
sider only H+Ba and He+Ba for G191-B2B (T e g = 60 000 K. 
log q — 7.6) and RE 0503-289 (7^ = 70 000 K, log g = 7.5), re- 
spectively. The Ba abundance was 1.5 x 10 -5 (mass fraction, 
about 1000 times solar, Asplund et al. 2009) in both mod- 
els. Figure 2 shows the respective Ba ionization fractions. 
Bavil is the dominating ionization stage in the line forming 
region (-2.5 < logm < 0.5). In the case of G191-B2B, 
the H+Ba spectrum shows two prominent lines in the FUV, 
Bavn /M943.10, 993.41A, and several much weaker Bavi and 
Bavil lines (Fig. 3). For RE 0503-289, Bavn A 993.41 A is 
the strongest line in the He+Ba, but some weaker Bavi and 
Ba VII lines are also visible. 

For the wavelengths higher than the FUSE wavelength 
range, an adequate observation is only available for G191-B2B 
(Sect. 3). Our test model shows a rich Ba VII spectrum within 
1150A < A < 1780A for this star, and a few weak BaV-Vl 
lines in addition (Fig. 4). The model spectrum for RE 0503-289 
shows the same lines, with deviations in the relative line 
strengths. We note that all these lines are much weaker than the 
four Bavi- VII lines in the FUSE wavelength range (Fig. 3). 

The preliminary models cannot be used for a reliable abun- 
dance determination because the neglected metals result in an 
unrealistic atmospheric structure (Fig. 5), and their missing 
background opacities have a strong impact on the strengths of 
the Ba lines. Therefore, we performed a precise determination 
of the Ba abundances based on detailed atomic data and elabo- 
rated model atmospheres (Sect. 6). 



1200 1300 1400 1500 1600 1700 

wavelength / X 


Fig. 4. Same as Fig. 3 {top: G191-B2B, bottom : RE 0503-289), to sim- 
ulate the HST/STIS spectrum (FWHM = 0.01 A). 


6. The photospheric Ba abundances in G191-B2B 
and RE 0503-289 

We determine the Ba abundances of G191-B2B and 
RE 0503-289 from their FUSE observations. Figure 6 shows 
that we can well reproduce the Ba lines with Ba mass fractions 
of 4.0 ± 0.5 x 10 -6 and 3.5 ± 0.5 x 10 -4 , respectively. These 
values are strongly over-solar (23 000 times and 265 times, 
respectively, Asplund et al. 2009) but in line with the determined 
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Fig. 5. Temperature structures of the final model atmospheres of 
G191-B2B {top panel, thick, red line; Rauch et al. 2013, H, He, C, 
N, O, Al, Si, P, S, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Zn, and Ge opaci- 
ties considered) and RE 0503-289 ( bottom panel, thick, red line; Rauch 
et al. 2014, He, C, N, O, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Zn, Ge, 
Kr, and Xe). The dashed, blue lines show the temperature structures of 
H+Ba {top) and He+Ba {bottom) models calculated with the same T e g 
and log g like the final models, respectively. The formation depths (i.e., 
r = 1) of the cores of our strategic of the cores of our strategic BaVI - 
VII lines in FUSE’S wavelength range are marked. 


abundances of other trans-iron elements (e.g., Werner et al. 
2012; Rauch et al. 2013, Fig. 7). It is worthwhile to note, that 
the Ba Vl/Ba VII ionization equilibrium is very well reproduced 
for RE 0503-289. 

In addition to the metals’ background opacity in the fi- 
nal model mentioned above, the determined Ba abundance of 
G191-B2B is only a fourth of that used in our test model 
(Fig. 4) and, thus, the lines in the wavelength interval 1 145 A < 
/l < 1750 A are weaker in the models, respectively. Figure 8 
shows a section of the spectrum where our test models pre- 
dicted the strongest lines (Fig. 4). A comparison to a model, 
that was calculated without Ba, shows that in this section only 
Bavn A 1472.96 A is a weak, isolated line that is visible in our 
model. Its strength is comparable to the low noise of the ob- 
served spectrum and, thus, we can determine an upper limit for 
the Ba abundance which is the same upper value determined 
from the FUSE spectrum (see above). For RE 0503-289, the de- 
termined Ba abundance is more than 20 times higher than in our 
test model (Fig. 4) and, thus, the potential identification of Ba 
lines in this wavelength range is most likely. A high-resolution, 
high-signal-to-noise HST/STIS UV spectrum of RE 0503-289 
is highly desirable. 

7. Results and conclusions 

We determined precisely the photospheric Ba abundances in the 
DA-type white dwarf G191-B2B (4.0 ± 0.5 x 10 -6 ) and in the 
DO-type white dwarf RE 0503-289 (3.5 ± 0.5 x 10 -4 ). These 
strongly supersolar abundance values are in line with those of 
other trans-iron elements in both stars (Fig. 7). 

Ba is one of the most massive s-process (slow neutron cap- 
ture) elements synthesized by low-mass (~1-3M 0 ) AGB stars 
(e.g., Lattanzio & Lugaro 2005). The s-process leads to abun- 
dance peaks at approximate mass numbers of 88, 138, and 208 



Fig. 6. Comparison of the theoretical line profiles of the strongest Ba VI 
and BaVII lines with the FUSE observation. The synthetic spectra were 
normalized to the observed local continuum fluxes. Top : G191-B2B, 
bottom: RE 0503-289. [X] denotes log ( mass fraction/solar mass frac- 
tion ) of species X. 
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Fig. 7. Photospheric abundances of G191-B2B (bullets, Rauch et al. 
2013, 2014) and RE0503-289 (squares, Werner et al. 2012; Rauch 
et al. 2012, 2014) compared with solar values (Asplund et al. 2009). 
The abundance uncertainties are about 0.2 dex in general. The arrows 
denote upper limits. The dashed line indicates the solar values. 


(West & Heger 2013; Karakas et al. 2014) due to closed neu- 
tron shells. Sr, Ba (isotopes from 128 Ba to 140 Ba), and Pb are 
located at these peaks and are used to represent the scaling of 
the s-process elements. Therefore, the Ba abundance determined 
in the hot WDs G191-B2B and RE 0503-289 establishes a new 
constraint for AGB and post- AGB stellar evolution and will help 
to understand the extremely strong metal enrichment. 

The identification of Ba lines in this paper was only possi- 
ble because reliable transition probabilities for Ba V, Ba VI, and 
BaVII were calculated. Analogous calculations for other highly 
ionized trans-iron elements is a pre-requisite for further iden- 
tifications an abundance analyses. The precise measurement of 
their spectra, i.e., their lines’ wavelengths and relative strengths, 
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Fig. 8. Section of the HST/STIS spectrum of G191-B2B compared 
with our synthetic spectrum. The thin (blue) model at the bottom is cal- 
culated without Ba. The identified Fe and Ni lines and the BaVII lines 
are marked at top. 

as well as the determination of level energies and the calcula- 
tion of transition probabilities remains a challenge for atomic 
and theoretical physicists. 
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